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ABSTRACT 

We present a critical review of theoretical models of Jovian 
decametric radiation, with particular emphasis on the lo-modulated 
emission. The problem is divided into three broad aspects: the 

mechanism coupling lo 1 s orbital motion to the inner exosphere, the 
consequent instability mechanism by which electromagnetic waves are 
amplified, and the subsequent propagation of the waves in the source 
region and the Jovian plasmasphere. At present there exists no comp- 
rehensive theory that treats all of these aspects quantitatively 
within a single framework. Acceleration of particles by plasma 
sheaths near Io appears to be a promising explanation for the 
coupling mechanism, while most of the properties of the emission may 
be explained in the context of cyclotron instability of a highly 
anisotropic distribution of streaming particles. The present state 
of the theory is evaluated, and some, suggested approaches for future 


work are discussed. 
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I. INTRODUCTION 

The nature and origin of the Jovian decametric radiation 
(DAM) have prompted much theoretical effort, especially since the 
remarkable modulating effect of lo -was demonstrated by Bigg [ 1964 ], 

The problem occupies an unusual niche in the gallery of radio- 
astronomical phenomena. As is true of emission from quasars, pulsars, 
and the sun, DAM presents a rich phenomenology. Unlike the cases of 
quasars and the sun, however, which exhibit detailed line spectra at 
optical wavelengths, radio observations [of both DAM and the weaker 
decimetric emission (DIM)] provided the only information about the 
plasma environment of Jupiter before the recent probes of the Jovian 
magnetosphere by Pioneers 10 and 11. Thus, for nearly twenty years 
the physical characteristics of the Jovian magnetosphere were inferred 
solely from radio astronomy. In this respect the problem is similar to 
that of pulsars. Like the sun, however, Jupiter is close enough that 
a detailed observational morphology of its radio emissions may be con- 
structed. This unusual combination of circumstances, together with 
the extensive development of understanding of terrestrial magneto- 
spheric processes, has allowed for wide variability in the mixture of 
inductive and deductive arguments in various individual attempts to 
explain DAM. 
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Despite the great amount of theoretical effort that has so far 
been expended , we have not yet arrived at a comprehensive, convincing 
understanding of DAM. Nevertheless, a review of the theoretical work 
to date is both pertinent and timely. The results from the Pioneer 
10 and 11 missions have excited widespread interest in Jovian physics 
in general, and much of the conceptual framework in which we evaluate 
these results has been developed in the course of attempts to eluci- 
date DAM. We anticipate that the relations between the problems of 
understanding DAM in particular and the Jovian magnetosphere in 
general will be more widely appreciated in the future. Furthermore, 
although we do not view any of the attempts heretofore to explain 
various facets of DAM as wholly success ful, we believe that in the 
paradigmatic nature of science a successful theory will emerge from 
the gradual modification and consolidation of previous work, rather 
than from radically new hypotl ses. In this spirit, vie attempt in 
this review to present a critical framework for the evaluation of the 
theoretical efforts to date, in the hope that such a framework will 
prove useful in both motivating and evaluating future work. The bulk 
of our discussion concerns the Io -modulated DAM; we discuss the lo- 
independent emission briefly near the end. 

A current ■review of DAM observations appears in this volume 
[Carr and Desch, 1975]. Previous reviews by Warwick [1964, 1967 , 
1970] and by Carr and Gulkis [ 1969 ] deal with observations of both 
DAM and DIM and the inferences about the Jovian magnetosphere drawn 
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from these observations. In addition, Warwick [1967] presents a 
critical review of some of the theoretical work up to 1967* 

Our approach in this review is to separate the theoretical 
problem into various more -or -less distinct facets, and to consider 
these facets both in themselves and in their structural relationships 
to each other. Most of the theoretical works to date deal principally 
with only one or two of the different elements of the problem, and so 
we discuss each work in the most appropriate contexts. We assume 
that the reader is familiar with the observations as they are 
described in the above-mentioned reviews, and so our references to 
observational results will usually be superficial and will seldom be 
referenced; the article by Carr and Desch in this volume will generally 
provide adequate observational background for our discussion. 

Some preliminary semantic definitions may be in order. As 
we have already done above, we shall sometimes distinguish between 
the phenomenology and morphology of DAM. By "phenomenology" we mean 
those aspects that are directly observed, such as intensity, time 
scales, bandwidth, polarization, and such features of the dynamic 
spectra as frequency drifts and modulation lanes. In contrast, we 
mean by "morphology" the organization of the data into such features 
as rotational profiles of occurrence probability and polarization, 
the similarity of the dynamic spectra of different events observed 
under the same geometrical conditions, the declination effect, the 
induction by Bigg of Io-modulation, the concomitant recognition of an 
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lo-independent component, and so forth. Moreover, we do not dis- 
tinguish between the familiar decametric emissions and the hectometric 
emissions observed from satellites by Brown [1974a, b] and by Desch 
and Carr U-974]. Except for the component of hectometric emissions 
that was designated by Brow [1974b ] as "mid -frequency" (MF) emission, 
the hectometric emission appears to be merely a continuation of the 
decametric spectrum to longer wavelengths, and we shall include the 
hectometric range in the traditional designation DAM. We note , 
however, that at this writing (mid- 1975 ) the observational statistic : 
are insufficient to determine whether either the "normal" or the MF 
components of hectometric emission are modulated by any of the Jovian 
satellites . 

Before plunging into the literature, we discuss some theoretical 
preliminaries . 

The high intensities, limited bandwidths, short time scales, 
and sporadic nature of individual DAM bursts indicate that they 
originate in stimulated emission from one or more collective micro- 
instabilities . The remarkable repeatability of the dynamic spectra 
observed under similar geometrical conditions, and the complex 
beaming pattern manifested in the source morphology, imply that the 
structure of the ambient medium is stable in the long-term sense and 
determines both the local conditions for plasma instability and the 
subsequent propagation of the radiation. Although this statement 
seems indisputable in the general sense, it must nevertheless be taken 
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with some qualification. From detailed study of certain repeatable 
features of dynamic spectra, Warwick [19^3] showed that these features 
always appeared for particular corresponding central meridian longi- 
tudes (CML), within a range of ± 9 °. Warwick interpreted these 
results as implying that the emission was beamed into a cone of 9 °. 
Alexander [1975], however, compared several pairs of identical- 
looking dynamic spectra observed one Jovian year (=* 11. 9 y) apart. 
Although the members of each pair of spectra commonly shared quite 
detailed features, Alexander concj.uded from a study of the rotation 
rates inferred from the different pairs that the beaming pattern 
might be subject to longitude shifts over a range of up to ± 10 °. 

Another point is so obvious that it is not often made explicitly, 
but we feel that perhaps for that very reason its implications are 
sometimes lost sight of in the literature. This point is that the 
occurrence probability distribution in the(CML, 7 Tq ) plane [ 7 ^ is 
the phase of lo from superior geocentric conjunction (SGC ) ] is con- 
ventionally interpreted as indicative of a beaming pattern. The 
occurrence probability over most of the plane, however, is substantially 
less than unity. This fact alone indicates that there must be 

some variable influence on the modulation mechanism or the waveguide 
structure of the plasma medium, or both. 

Plasma instabilities require a source of free energy. Large- 
amplitude magnetohydrodynamic (MED) waves or shocks may provide free 
energy under some circumstances. In micro-instabilities (by which 


8 


we mean instabilities that are predicted only in kinetic-theoretical 
analyses and cannot be predicted from fluid theory) the free energy 
resides in some nonthermal feature of the particle distribution 
function F(r, v, t). Such a nonthermal feature may be manifested in 
either the configuration-space or velocity- space dependence of F . 

For -ixample, density gradients in a magnetised plasma create drift 
currents , which may drive a variety of instabilities. Velocity-space 
features such as beams of high-velocity particles, gaps in the 
distribution in some velocity range, or anisotropy in either pitch 
angles or Larmor phases, may result in instability. 

The repeatability of DAM spectra suggests that the instability 
occurs for waves of frequency a,t or related to one of the character- 
istic frequencies of the plasma, as is commonly the case. In the 
high-frequency regime, these are the electron plasma frequency f , 
the electron gyrofrequency f , and the upper hybrid frequency 
f UH = ( f pe + f c e)"^ 2 ' ^ Jupiter, the DAM frequency range may be 

identified with these haract eristic frequencies only near the planet: 
roughly speaking, at altitudes of ^ 1 Rj. 

Thus , the central problem in the theory of DAM is to elucidate 
the coupling mechanism whereby the orbital motion of lo provides a 
free-energy source in the inner plasmasphere . The existence and 
nature of an instability are determine . by the form in which the free 
energy appears . The coupling mechanism and the subsequent instability 
mechanism by which the radio emission is amplified are, however, viewed 
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here as two logically distinct elements of the problem. Roughly 
speaking , we might expect that the coupling mechanism primarily 
affects the DAM morphology, while the nature of the instability is 
reflected in the phenomenology. This statement is intended simply 
as a general guide, however, and not as an organizing postulate. 

The third element of the problem is the propagation of the 
radiation in the Jovian plasmasphere. Propagation effects are 
expected to be manifested in both the phenomenology and morphology. 

As we shall see below, for example, certain instability mechanisms that 
have been proposed actually amplify waves that are trapped in the 
Jovian plasmasphere. Therefore, some mode-conversion process is 
necessary in order for the radiation to propagate into free space. 

In this case the polarization of the escaping wave, which is generally 
the datum by which observers identify the base mode, may be determined 
by the escape mechanism rather than by the instability. The escape 
mechanism may also determine the initial wave-vector of the escaping 
ray, thus deter mi ning the ray path in the plasma medium and the 
beaming pattern. 

Provided the boundaries between the different aspects of the 
problem are not regarded as impenetrable, our formulation is a useful 
one which enables us to construct a general conceptual framework 
within which to consider particular theories. 
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II. THE COUPLING MECHANISM 

The question of the coupling mechanism has "been addressed, 
to varying degrees, hy many authors. Smith and Wu [1974] divided 
the majority of the proposed mechanisms into three broad categories: 

(i) Wave models , In which the motion of a conducting Io 
through the magnetoplasma may generate whistlers or MHD waves that 
propagate to the Jovian ionosphere; 

(ii) Acceleration models , in which the electric field 
induced in the frame of Io "by the satellite's motion through the 
ambient magnetic field is invoked to accelerate particles to 
energies from tens to hundreds of keV. These particles then 
precipitate along the lo-threaded flux tube (IFT); 

(iii) The sweeping mode l, in which energetic particles 
trapped on the IFT are removed from the trapping volume by 
impingement on Io. 

To our knowledge, the above categories embrace all the 
cases in the literature in which a coupling mechanism has been 
examined quantitatively. This does not mean that their corresponding 
instability mechanisms have been considered quantitatively; as we 
shall see in section III, no quantitative model of an instability 
mechanism that migno result from MHD-wave coupling has yet been 
presented in the literature. Whether or not such mechanisms may 
exist, the generation of MHD-wave disturbances by Io may have 
important consequences which are discussed below. 
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In addition to the coupling models described above., other 

more speculative hypotheses have been advanced. Duncan [1970] 

suggested that the source morphology was determined by a combination 

of the geocentric beaming pattern and heliocentric effects. He 

constructed a model of the Jovian magnetosphere in which the 

subsolar point of the magnetopause was at 5 Rj, and conjectured 

that the strong confinement of Io modulation to orbital phases 

near 7 ^ — 90° SGC and 240° SGC implied that Io crossed the 

magnetopause at about these phases. Duncan further speculated 

that these magnetopause crossings would precipitate trapped 

particles into the Jovian ionosphere, although he did not indicate 

why or how this might occur. To evaluate Duncan's hypothesis, however. 

it suffices to note that the Pioneer 10 and 11 probes have revealed 

the Jovian magnetosphere to be extended to characteristic 

dimensions of 50-100 R in the dayside hemisphere. 

d 

Another phenomenological suggestion in a spirit si mi lar to 
that of Duncan's was made by Conseil et_ al. [1971/ 1972]. These 
authors inferred a correlation between the velocity of the solar 
wind at Jupiter (extrapolated from satellite measurements at 1 AU) 
and the phase of Io at the onset of DAM storms. They suggested 
that the correlation might be attributed to a plasmapause bulge, 
similar to that of the terrestrial plasmapause, the location of 
which varied in response to solar wind ram pressure. They 
speculated that Io might trigger DAM storms upon passing through 
the plasmapause boundary. 
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The statistical significance of the correlation asserted 
by Conneil eb al. is unclear; it has been questioned by Carr and 
Desch [1975]* Furthermore, a Jovian plasmapause would be 
unlikely to contain a bulge analogous to that in the terrestrial 
case. It is believed that the terrestrial plasmapause occurs 
at the boundary between closed and open streamlines arising from 
the combined flow fields of corotation and convection. 

In the Jovian magnetosphere, corotation would be expected to dominate 
convection out to large distances, because of the rapid rotation rate of 
Jupiter. A comparison of the flow pattern to be expected at Jupiter 
and that in the terrestrial magnetosphere is given by Brice and 
loannidis [1970]. Although this picture may be subject to 
considerable modification in light of the unexpectedly large 
extension of the Jovian magnetosphere and the dominance of its 
outer structure by plasma outflow, it should not be greatly changed 
in the inner region. A true plasmapause at Jupiter would probably 
be determined by the breakdown of corotation. 

Frank et. ajL [1975], however, inferred the existence of a 
plasmapause at the orbit of lo, where the density of protons of 
energies 100 eV - 4.80 keV was observed to increase by a factor 
of about 3 to 4 on the inbound pass of Pioneer 10. As may be 
seen in Figure 1 [Frank et al« , 1975] > however, this density 
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decreases again inside Io’s orbit, before resuming a general increas- 
ing trend as the radial distance decreases. These observations 
could also be explained by injection of these particles by Io, rather 
than as a plasmapau.se structure. 

Whether or not the structure of the thermal-plasma density 
near Io is indicative of a plasmapause, it is clearly complicated. 

The basic idea behind the suggestion of Conseil et al - was that the 
stimulation of DAM by Io might depend on the local plasma density. 
This density may, in fact, be an important parameter in the MHD-wave 
model of Goertz and Delft [1973; c.f. also Deift and Goerts, 1973; 
Goertz, 1973a, b] and the recent acceleration models of Shawhan 
et al. [19743 and of Smith and Goertz [1975 3. 

We now consider the quantitative categories (i) - (iii). As 
was mentioned above, in most of the works discussed in this review 
the major focus is on only one aspect of the total problem; in the 
majority of the papers we consider in this section it is the coupling 
mechanism. In most of these papers, however, the authors had some 
hypothesis in mind concerning the subsequent generation or propaga- 
tion of the electromagnetic waves. In some other works, such as 
that of Wu [19733, two or more aspects of the problem are treated 
integrally. Therefore , in order not to fragment our exposition 
unduly, we shall occasionally anticipate some points of the later 
discussion in order to place the immediate topic in context. 
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A. Wave Models 

Ellis [1965] first speculated, that the coupling might he 
due to either MHD or "whistler ("which he called "electromagnetic”) 
-waves generated by lo. He did not investigate how such waves 
might he excited, except to note that the relative velocity of 
Io normal to the corotating magnetic field was greater than the 
phase velocities of these modes, all of which have resonances 
(wave number k •*> <=) at angles 0 < rr /2 with respect to the field, 
Ellis, therefore, assumed that these low-frequency modes might 
he generated hy Cerenkov radiation from Io. This is by no means 
obvious, however, because in order to he a source for the waves, 

Io must either have an appropriate current distribution on its 
surface or cause such a current distribution in the ambient plasma 
by perturbing its flow. 

Ellis hypothesised that the Io-generated waves would 
interact with pre-existing particle streams or "bunches” to 
stimulate the electromagnetic radiation observed as DAM, in 
analogy to the well-known phenomenon of stimulated VLF emissions 
in the terrestrial magnetosphere. Thus, the whistler waves are 
the important ones in this context; this mechanism was also 
considered by Chang [1963], before Io-modulation had been recognized. 
It is implicit in this scheme that the resulting DAM waves travel 
along field lines, and that the source region would therefore 


be the CML. 
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McCulloch [ 1967 ] developed Ellis’s ideas for the coupling 
more quantitatively. He too assumed whistlers to he generated 
hy Io, and did ray tracing in a model magnetosphere. By varying 
the parameters of the model, and hy assuming that the coupling 
region where the whistlers interacted with particle streams was at 
magnetic latitudes between 75° and 80°, McCulloch obtained a 
distribution of coupling regions in the (CML, 7 j q ) plane, for 
which the relative intensities of the propagated whistlers 
corresponded well with the distribution of occurrence probability 
of BAM (c.f. in particular his Figures J, 9? amd 10). In 
particular, McCulloch’s results indicated that the time required 
for the propagation of the whistlers from Io to the Jovian 
ionosphere was of the order of 55 min. This is the time required 
for Io to move 15° in orbital phase relative to the corotating 
medium; in this way, McCulloch accounted for the asymmetry of the 
Io phase, for peak occurrence probability, about the Earth- Jupiter 
line (Figure 2). 

The parameters which McCulloch found necessary to describe his 
model magnetosphere, in order to obtain good agreement between the 
coupling mechanism and observational morphology, are not in accord 
with the Pioneer observations. For example, McCulloch used an 
ionospheric electron temperature of 2000 K, whereas KLiore at a.1 . 
[19753 have inferred T == 750 K. Furthermore, McCulloch states 
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[1967, p* 1502] that he uses a density model for the thermal plasma 

of the form derived by Melrose [19671- But he also states that "(the) 

most satisfactory models, however, were those for which the electrons 

were concentrated in a disc about the magnetic equator, and the plasma 

frequency and gyrofrequency at lo -were nearly equal" [McCulloch, 

1967, p- 1305 3- The latter condition, f cg =- f , is actually in 

accord with observations: for a dipole moment of 4 G-Rj^ [smith. ~et al. , 

19753, the field strength B q = 0.02 G at Io and f - 56 kHz, while 

for a plasma density at Io of =* 30 cm - ^ [Frank at al. , 19753, 

f =* 55 KHz- The feature of a plasma disk in the equatorial plane, 

pe 

however, is not contained in Melrose's density model, in which the 
density decreases monotonically with latitude along a field line. 
Therefore, the possibility exists that McCulloch’s work may be inter- 
nally inconsistent. Finally, McCulloch’s ray tracing calculations 
assumed a centered, tilted dipole field of equatorial field strength 
15 G ; it is not clear in context (c.f. p. 1306) whether ad hoc 
anomolies, such as were postulated by Ellis and McCulloch [ 1963 3 
were incorporated or not. As has been shown by Acuna and Ness [19753 
and Smith et al. [ 1975 3 > however, the field near the surface is 
strongly influenced by higher-order moments. 

■Hr 

The plasma disk is contained in density models developed by 
Gledhill [19673 and ‘Goertz [1973a, 19751- 
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Apart from questions about the associated instability mechanism 
and consequent propagation effect s, these considerations indicate that, 
at the very least , coupling by -whistlers must be re-examined in the 
context of the plasma environment deduced from the Pioneer observa- 
tions. This comment does not apply uniquely to the work of McCulloch, 
however . 

Other authors have concentrated on low-frequency (f « f ) MED 
waves to couple lo to the DAM source regions. Marshall and Libby [ 1967 ] 
suggested that such waves might induce spin transitions in free radi- 
cals in the Jovian ionosphere, presumably because of the presence of 
off-diagonal terms in the pressure tensor, in analogy to the inter- 
action of phonons with free radicals in crystals [Warwick, 1967]. 

The first quantitative estimates of MHD-wave coupling were made 

by Warwick [ 1967 ]* who considered lo to be perfectly conducting and, 

therefore, argued that it would push aside the ambient field. He 

estimated that for "resonant 1 ' waves of phase velocity V ^ = V^ o 

(where = 56 km/sec is the speed of lo relative to the corotating 

plasma) and wavelength 2 R^, waves of amplitude B = B q 

might be produced, where is the Alfven velocity near lo. At the 

current estimate of V. == 700 km/sec, this would give B U B . From 

other considerations, involving large but finite conductivity for lo, 

Warwick [1970] estimated B =- 2.6 B . Waves of this magnitude are 

wo “ 

inherently nonlinear, and would be expected to steepen into shocks 
as they propagated toward Jupiter. Warwick [1970] estimated that 
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about 4 X 10*^ watts would thus be radiated from Io. He did not 

attempt to partition this energy among the various MHD modes ( shear 

Alfven, congressional Alfv£n, and magnetosonie) . Instead, he noted 

o 

that if the MHD radiation from lo were isotropic, about 5 X ICr W 

c o f 

would be intercepted by an area of 1.6 X Hr km. (the presumed size 
of the DAM source region) at the surface of Jupiter. 

The estimates by Warwick [ 1970 ] were made partly in the context 
of a critique of the paper by Goldreich and Lynden-Bell [19693* which 
is discussed below. These latter authors assumed that the electric 
field induced across Io by its motion relative to the Jovian magnetic 
field would drive a current system that closed ir the planetary iono- 
sphere. In contrast, Warwick asserted that the current system would 
be closed within Io itself. This idea was also adopted by Schmahl 
[19703? "who considered Io to have a permanent dipole moment in- 
duced by the DC average of* its ambient magnetic field, and thus to 
be a point current source 

D s (r, t) = c M IO X V 6 (r - V Iq t) 

Schmahl calculated that 8 X Hr W would he delivered to the Jovian 
ionosphere in the shear-Alfven mode, which propagates essentially 
along the IFT. 

In a contrasting approach, Goertz and Deift [19733 assumed 
that currents leak from Io into the ambient plasma. They attempted to 
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calculate the steady flow pattern of the ambient medium around Io. 

The resulting distortions of the ambient magnetic field are to be 
understood in the sense of MHD wave fields. Some particularly strik- 
ing features of their model are: (i) the perturbations of the field 

ahead of Io propagate to the Jovian ionosphere as solitons and are 
reflected as solitons, leading to a stretching of the field roughly 
analogous to the stretching of a string; (ii) the existence of an x- 
type neutral point behind Io (Figure 3)* The x-type configuration is 
asserted to occur when the field lines leading Io are stretched too 
far; reconnection then takes place behind the moon, leading to the 
propagation of Alfven waves up the IFT. In effect, Io "plucks” the 
magnetic field lines, like plucking a violin string. The angular 
limit of distortion before reconnection occurs is about 15 °, account- 
ing for the lo-phase asymmetry. Goertz and Delft estimate the maxi- 

* 15 

mum power radiated in Alfven waves to he s 2.7 X 10 ' W. 

In Goertz and Delft's model, the wavelength of the Alfvdn- 
waves generated hy reconnection is taken to be 6 Rj 0 (i.e. , k= R~ q ). 
The resulting frequency of the waves is between 1-3 Hz, depending on 
the plasma density at Io. This density is determined from a magneto- 
spheric model developed by Goertz [1973a]. Delft and Goertz [ 19731 
consider the propagation of these waves along the IFT. They find a 
"transmission coefficient" for the waves into the Jovian ionosphere: 
the depth to which the Alfven waves penetrate— and, therefore, the 
gyrofrequency at this depth — is a function of the wave frequency. 
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Goertz [19755.1 equated the gyrofrequency at the local altitude of 
penetration, as a function of longitude, to the maximum DAM frequency 
at that longitude, showing good agreement with observation [Goertz 
and Has chick, 1972 J- 

Goertz and Deiffc attempted to consider the three-dimensional 
nature of the MHD flow around Io. It is not clear, however, that 
they did so consistently. They assumed the conductivity of the ambient 
medium to be a scalar, whereas In a low-P plasma (p == OrrM^/B ) 
the tensor properties of the conductivity are significant. Further- 
more, even within the context of their model, the field lines which 
obey the geometry of Figure 3 are limited to those lines near the 
0 stagnation point of the flow; in principle, these lines degenerate 
to a set of measure zero, and it is then unclear to what extent their 
analyzed flew Is dominated by the actual three-dimensional flow. In 
particular, the efficiency of flux reconnection at x-type neutral 
points may be drastically lower in three-dimensional flow than in two- 
dimensional flow. Furthermore, the back- flow towards Io associated 

■X- 

with the reconnection is not taken into account. 


In defense of these authors, it must be noted that the three- 
dimensional MHD flow past a spherical body is an intractable problem, 
particularly in a low-(3 plasma, nevertheless, the spirit of our 
comments is that a two-dimensional solution may be greatly misleading. 
In the common example, a tight -rope walker confined to the vertical 
plane could never fall off the tight -rope. 
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In addition, Deift and Goertz consider incompressible plane 
waves propagating along the IFT. It is not clear whether, in the 
situation of wave generation by flux reconnection, compressible 
modes or finite wave packets [Goldstein, Klimas, and Barish, 1974] 
might be more appropriate. 

Finally, we remark that the entire discussion by Goertz and 
Deift [19753, Deift and Goertz [1975], and Goertz [1973b] is heavily 
dependent for its quantitative significance upon the magnet ospheric 
model of Goertz [1973a]. 

Hone of the wave-coupling models discussed include a 
quantitative discussion of a consequent instability mechanism. In 
respect to the instability mechanism, however, the power available 
in the coupling mechanism is of vital importance because the 
efficiency of conversion to electromagnetic waves is necessarily low, 

g 

We shall take 10 W to be a nomimal power requirement for electro- 
magnetic waves in DAM bursts [Warwick, 1970 - 1 . Therefore, the power 
in the coupling mechanism must be » 10® W. 

McCulloch [ 1967 ] estimates that 2 x 10^ W may be transferred 
from lo via whistlers. Taken by itself this figure would indicate 
a conversion efficiency of 0. 05, which is incredibly high. The 
presumed interaction mechanism, however, is the phase-organization 
of trapped particles (c.f. section III), and therefore the particles 
may be more important in the energy balance than the whistlers. 
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Because no instability mechanism has been demonstrated for 

coupling by low-frequency MHD-waves, we take the viewpoint that 

the conversion efficiency must be compatible with the power 

available in the MHD waves. (This viewpoint is consistent with 

the suggestion by Goertz [1975b] that the DAM radiation may 

originate from particles accelerated in the MHD-wavefronts . ) In 

10 .. 

this light , Schmahl T s estimate of = 10 W in shear-AlfVen waves 

_2 

requires a very high conversion efficie. cy of 10 . On the other 

15 

hand, the maximum-power estimates of 2.7 X 10 W by Goertz and 
l6 

Deiffc and 4 X 10 W by Warwick [1970 1 are much too high to agree 
with the recent determination by Goldstein [1975 ] that the secular 
acceleration of Io is limited to a value 

S < 3 x 10" 18 sec” 1 , 

P 

where p is the orbital period of Io. Because angular momentum 

is conserved in the system of Jupiter and its satellites, this 

14- 

finding limits the power dissipated by Io to about 10 W. Although 

such a power is consistent with conversion efficiencies for DAM 

-6 

bursts of £ 10 , it also indicates that the coupling mechanisms 

considered by Warwick [1970] and by Goertz and Deiffc c ann ot work 
at the maximum powers estimated by their authors. 

It is certain, however, that Io must create some MHD 
disturbance in the magnetospheric plasma. Whether or not these 
waves are responsible for the coupling to DAM, they may have other 
significant effects. One such effect could be to provide local 
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ionospheric heating near the base of the IFT, A crude calculation 

(Appendix A) yields an estimate that the waves might heat the 

4 5 

ionosphere to temperatures of the order of 10 to 10 K. Because 
the plasma temperature may, under some conditions , be an important 
parameter in the contexts of instability mechanisms and wave 
propagation, such heating may have significant consequences for 
these aspects of the theory. 

B. Acceleration Models 

The starting point for acceleration models is that, owing 

to the motion of Io with velocity jv^| = 56 km/sec relative to the 

-+ *+ -> , 

corotating plasma, an induced electric fie] 1 E in(i = Vj X B /c 

■+ 

appears in Io's frame, where B^. is the ambient magnetic field at Io. 
If Io is not a good conductor, then this field would he attenuated 
only by volume polarization throughout its interior. Piddington 
and Drake [ 1968 ], however, suggested that Io might be highly 
conducting and so its interior would be screened from the induced 
field by surface polarization charges. As a result, they asserted 
that the induced emf of 400 kV could be mapped onto the field lines 
of the IFT, which are assumed to be equipotentials. The existence 
of the potential is communicated along the field line at the Alfven 
speed = B q /( 4 ttp^) s , where B 0 Is the magnetic field strength and 


p^the mass density of the plasma. Provided that the conductivity cr 

is large enough so that the field does not diffuse through Io in a 

time shorter than the propagation time of Alfven waves to the Jovian 

ionosphere, a DC current circuit which closes in the Jovian ionosphere 

will he established. The resistive diffusion time j , over a scale 

nr 

length by 


m 


- h-TT 


E 


Io 


The precise value of is not crucial in the present context. For 
purposes of estimating a constraint on a, we may set =“ 60 sec; this 
is the time required for Io to move a distance equal to its own dia- 
meter. Then for > 60 sec, we require <j > 1.5 X 10 7 sec" = 1.7 X 30 _? 
mho/m. 

Goldreich and Lynden-Bell assumed sufficient conductivity for 

Io that the DC circuit model would be valid. They asserted that this 

requires T > 2r fl as we have defined r A . In actuality, the DC circuit 
m a A 

will he established for ~ r^, because only a one-way trip for an 
Alfven wave is necessary to establish the potential In the ionosphere. 

Strictly speaking, it is communicated at the phase velocity m/k 

v- p 2 ^ 

of an Alfven wave, where m/k = V ^ = c/( 1 + c/V^) 2 . For « c, 

Vpk ■+ V^; in the limit V ^ -* c. The opposite limit, Y^ 0, 

cannot he recovered from the Alfven-wave dispersion relation. 
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(One may turn on a light switch and get light very quickly; no electron 
in the wires need traverse the entire circuit.) Goldreieh and Lynden- 
Bell, however, actually required more than the existence of a circuit 
which closes in the Jovian ionosphere. They argued that if r » 2 
the entire IFT will be frozen to lo, only slipping in the ionosphere 
in order to follow rigidly Io's orbital motion. The importance of 
this assumption to Goldreieh and Lynden-Bell was that it enabled them 
to argue that the field vanishes within the flux tube and that 
the current is carried in thin sheets on the boundary of the IFT. 

From this argument, and from some assumptions regarding the thickness 
of these sheets and the plasma density, they inferred a requisite 
energy for the current-carrying particles of a few keV. 

Goldreieh and Lynden-Bell assumed that there would be some 
potential drop along each field line. They then solved the 
current circuit equations to find the shape of the potential 
distribution across the IFT. Their solution has the correct form 
if one assumes that the current flow may he proportional to the 
voltage [Smith and Goertz, 1975]- For the potential to drop along 
the field line, however, requires some parallel resistivity. 

Goldreieh and Lynden-Bell did not discuss the origin of this 
resistivity, and so their solution for the potential was left in 
terms of an unspecified constant V . 

The models by Piddington and Drake and by Goldreieh and 
Lynden-Bell are oversimplified in their assumptions regarding 
the MHD interaction of Io with the ambient plasma. Gurnett [1972] 
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considered, the effect of the formation of potential sheaths 

around. Io. An object placed in a plasma attains a surface charge 

because the flux of thermal electrons hitting its surface is 

greater than the flux of iors. As a result, there develops a 

negative surface potential •which is shielded from the ambient 

plasma by a Debye sheath. The equilibrium state is one in which 

the current to the surface is zero. If the surface emitB a 

photoelectron current larger than the current it collects from the 

thermal flux, the potential at the surface becomes positive and 

a photoelectron sheath of the opposite polarity develops. Gurnett 

showed that in the ease of a sunlit Io, part of which emits 

photoelectrons and part of which is in darkness, both types of 

sheaths occur. Their respective current densities -T^ (for the 

Debye sheath) and J (for the photoelectron sheath) and areas 

P 

and are related by the current balance condition: 


JA = 
P P 


J D A D 


Gurnett ' s model was considered in more detail by Hubbard, 
Shawhan, and Joyce [1974], and was extended to account for the 
existence of an ionosphere at Io by Shawhan et al. [1974]. The 
principal features of the latter modification are: (i) the conduc- 

tivity required to close the current system is provided by the 
ionosphere of Io, rather than by the solid moon; (ii) the sheaths 
are assumed to form at some ill-defined height at the "top" of 
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the ionosplK r The electron density profile of lo's ionosphere 
deduced hy KLiore et al. [1974] shows a precipitous drop at an alti- 
tude of about 750 km, where 'che scale height changes abruptly from 
approximately 220 km to approximately 36 km. Shawhan et al. inter- 
preted this as the sheath region. Adding this altitude to the radius 
of lo to give an "effective" radius for the conducting system, the 
total induced emf becomes 370 kV. The sheath model also assumes a 
current system that closes in the Jovian ionosphere j the model Is 
schematically depicted in Figure 4. 

Gurnett [1972] considered an idealized steady-state model of 
the current system (Figure 5 ) and concluded that If the height inte- 
grated Pedersen conductivity ^ of the Jovian Ionosphere is lower 
than = 5 mho, the I FT would move with Io and the bulk of the potential 
drop in the external circuit would occur In the Jovian Ionosphere, 
rather than in the sheaths. Conversely, Gurnett asserted that if 

E ^ 3 mho, the IFT slips by Io and substantial voltage drops occur 
P 

across the sheaths. After noting these conclusions, Hubbard et al. 
[1974] assumed that the major potential drops occurred across the 
sheaths [c. f. their Eq. (6)3. This assumption was also made by 
Shawhan et_ al. [1974] and by Shawhan [19751* All of these authors 
assumed that particles traversing the sheaths would gain potential 
energy equal to the sheath potential. 

Gurnett estimated J^ and J^ independently, as did Hubbard 
et al. Similarly, Hubbard et_ al. estimated the sheath thicknesses as 
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independent parameters , and similar estimates were made by Shawhan 

et al. These latter authors, however, recognized the sheaths to he 

potential double layers, which have been investigated by Block [1972], 

Carlqvist [1972], Knorr and Goertz [197*5-], and Goertz and Joyce [1975]. 

Such layers form when electrons and ions in a plasma with T g ^ T. 

drift relative to one another with a velocity U ^ V g , where V g is the 

electron thermal velocity. A potential then develops over a scale 

length L. Smith and Goertz [1975] noted that the development of the 

double layer was subject to the Buneman instability [Buneman, 1958], 

which produces acoustic turbulence in the frequency range ^ « a, « m . 

Tne convection of the turbulence at the group velocity, which is 

given by V =< 0.4 U for the most unstable waves, leads to a natural 
gj-D 

scale length L for the sheaths, given by 


3 

B 


where y =* ® Q /l8 is the maximum growth rate of the Buneman instability. 

Smith and Goertz also showed that the limiting drift velocity U was 

related to the sheath potential § by 

s 


1 

2 


m 


u 2 



By considering a steady-state current system similar to that analyzed 
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by Gurnet t, Smith and Goertz showed that the maximum potential drop 
across the sheaths is limited to about 140 kV. 

We stress that all of these analyses regarding the sheath 
model assume the current circuit to be closed in the Jovian iono- 
sphere. The self-consistent system of drift currents in the iono- 
sphere of Io, the sheath regions, and the ambient plasma into which 
Io is penetrating (and which is therefore magnetodydrodynamically 
perturbed) is not analysed, and it has not been demonstrated con- 
clusively that the current system does not close locally near Io. 

We note, however, that the current system cannot close in a region 

— > 

which moves with Io, because in such a region V x j = 0. 

Related to this point is the question of whether the conduc- 
tivity of Io’s ionosphere is itself sufficient to carry the current 
without causing a large drop of the emf. Shawhan et al. [1974] 
estimated the conductance of Io's ionosphere to be 260 mho, a value 
sufficient by more than an order of magnitude. The atmospheric model 
of McElroy and Yung [1975] is also compatible with a conductance of 
this order. This model assumes and ML to be the primary atmo- 
spheric constituents. Whitten, Reynolds, and MIchelson [1975] 
developed a model in which the primary constituent is neon. In this 
model, the ionospheric structure inferred by KLiore et al. [1974] 

Is obtained with a much smaller neutral density than in McElroy and 
Yung’s model, and the ionospheric conductivity is correspondingly 
lowr, perhaps becoming marginal with regard to the requirements of 
the sheath model [Shawhan, 1975]- 
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C. The Sweeping Model 

The MHD-wave and acceleration models discussed above all entail 
the assumption that Io is a good conductor. In contrast, Wu [1973] 
assumed that the conductivity of Io was sufficiently low that the 
Jovian magnetic field would diffuse readily through it. In that 
case, particles trapped on the field lines threading Io would impinge 
on the moon and he "swept” from the magnetosphere. Sweeping had 
been previously discussed by Mead and Hess [1973] in the context of 
its effects on the structure of the Jovian radiation belts. 

Energetic particles trapped in a dipole magnetic field execute 
a complete bounce, from a given mirror point to the conjugate mirror 
point and back again, in the time 

R. 

t. =* (l.3 - O .56 sin a ) 4L — , 
b o v 

where is the equatorial pitch angle, L is the conventional mag- 
netic shell parameter, and v Is the speed of the particle. For 
particles at L a* 6, mirroring at gyro frequencies corresponding to 
the DAM frequency range, ci^ ^ 0.1 rad, and so ^ 30 Rj/v. In the 
frame of reference corotating with Jupiter, Io moves a distance 

equal to its own diameter in a time t — 60 sec. Therefore, particles 

s 

that execute a half -bounce in a time shorter than t will impinge 

s 

upon Io and be removed from the field lines described by their 
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guiding center motions. The threshold velocity v c above which 
particles will he swept in this way is determined hy equating 
T b (v )/2 = t . For L = 6, =* 1.8 X 10^ km/sec, corresponding to 

threshold energies of approximately 0.8 keV for electrons and 1.6 
MeV for protons. Particles above the appropriate threshold are 
effectively swept from the lo flux tube. 

In contrast to the instability mechanisms associated with 
other coupling mechanisms, that due to the sweeping effect depends 
on the energetic protons, rather than electrons. In section III we 
shall consider the form of the proton distribution function; here 
It suffices to note that it includes as a parameter a characteristic 
energy E . 

ir 

The sweeping of energetic particles by lo results in a cavity 
in the density distribution of such particles. In the orbital plane 
of lo, the cavity has a cross-section similar to a wake and remains 
for some time after the passage of lo; diffusion across the cavity 
boundary tends to replace the energetic particles on some time scale 
T-j-j, long compared to the growth time of the electromagnetic insta- 
bility. The cavity cross-section in Io T s orbital plane is illustrated 
in Figure 6. The density of. particles with energies above the 
sweeping threshold falls to zero across the distance 6. The bow 
portion of the cavity where the sweeping is incomplete has a char- 
acteristic dimension f given hy 
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& = 




V, 


lo 


thus the how of the completely swept portion has a cross-section 
which is a semi-circle displaced from the how of lo by the distance 
j L The dimension 5 of the lateral boundary layer is given by the 
intersection of this profile with the trailing limb of lo: 


8 = IL. 
To 



cos sm 


-1 



(II-1) 


Equation (IX-1) is valid as long as it gives a value for 6 greater 

than the characteristic proton Larmor radius R = (E /MQ ) 3 , where 

P P‘ P 

fip = 2rrf c p is the angular proton gyrofrequency; B is the physical 
lower bound on 8. The equatorial boundary thickness 8 maps to values 
6 '(A.) < 6 at magnetic latitudes |l| >0 (Figure 7): 


s '(>0 “ 8 Jr ®.' • 

Assuming R (o)/s « 1, we have R (>-)/&' < 1 also, and the boundary 
P P 

thickness at higher latitudes is given approximately by the mapping 
for a centered collinear dipole: 


5' _ cos^ \ 
6 \k - 3 cos 2 
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Owing to the density gradient of energetic protons across the 
walls of the IFT, drift currents appear. These currents drive insta- 
bilities at the upper and lower hybrid resonances; these instabilities 
are discussed in later sections. 

D. Evidence From Pioneers 10 and 11 
The trajectories of Pioneers 10 and 11 did not pass through 
the instantaneous IFT, and so the evidence from these probes regarding 
the coupling mechanisms discussed above is indirect and inferential. 
Inbound, Pioneer 11 passed about 6000 km from the instantaneous 
IFT as defined by the D^ field model of Smith et al. [1975 1 at 
about 0300^ on 3 December 197^- [Fillius, Mcllwain, and Mogro-Campero 
1975]. The spacecraft was then at approximately -40° magnetic 
latitude, and so its distance from the IFT was about twenty times 
the radius of the flux tube itself at that latitude. No distortions 
of the magnetic field near this region have been reported in the 
literature, indicating that the MHD perturbations arising from 
Io T s penetration of the magnetospheric plasma are either unre solvable 
or are localized at this latitude to distances much less than 
6000 lan from the IFT. Thus, it is unlikely that To warps its 
flux tube by 15°, as suggested by Goertz and Deift. No direct 
evidence about the actual MHD nature of the lo-magnetosphere 
interaction can be obtained from trajectories such as those of 
the Pioneers, however. 
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The particle data from both Pioneers 10 and 11 show that 
sweeping occurs and is effective. In this regard, it must 
be remembered that sweeping of the IFT is followed by diffusion 
of the energetic particles back into the swept region. The 
Pioneer 10 data were often interpreted to indicate that sweeping 
is much less effective than predicted by Mead and Hess] typical 
reductions In the particle fluxes at the IFT observed on Pioneer 
10 were by factors of about two to five, and the sweeping effect 
seemed generally to be confined to lower-energy particles. Mead 
and Hess, however, considered only the slow processes of electric 
and magnetic diffusion, as they are understood to operate in 
the terrestrial magnetosphere, to balance the sweeping effects 
of the various Jovian satellites. They were led to predict 
large decreases, by several orders of magnitude, of the trapped 
particle fluxes at the orbits of each of the moons. Smith [1973], 
however, showed that the sweeping effect leads not only to the 
high-frequency electromagnetic instability described by Wu, 
but also to a low-frequency instability which produces strong 
diffusion ocross the boundaries of the IFT. (To the extent that 
sweeping is effective at other satellites, the same mechanisms 
operate there also. ) Moreover, the diffusion coefficient is 
energy -dependent, increasing with the particle energy. Thus, 
the relative mildness of the sweeping effect as observed in the 
Pioneer 10 data may be attributed to the fact that the satellite 
crossed the orbit of Ic several hours behind the satellite, and 
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the density gradient had been eliminated in the case of higher- 
energy particles and considerably reduced at lower energies. 

Further evidence of this interpretation may be found in the Pioneer 
11 data. The trajectory of Pioneer 11 was such that the spacecraft 
passed the lo orbit much closer to Io, both inbound and outbound, 
than did Pioneer 10; sweeping was found by all the particle experi- 
ments to be quite severe. 

There is clear evidence for acceleration of electrons at 
Io, but its interpretation in terms of existing models is ambiguous. 
On both the inbound and outbound crossings of Pioneer 10, Fillius 
and Mcllwain Ll9/4l observed enhancements, by a factor of about 
two, in the electron flux at energies E ■ > 0.16 Mev. These enhance- 
ments were observed just inside the orbit of Io, in accord with the 
sheath-model predictions of Shawhan et. al L 19 733. On the inbound 
crossing of Pioneer 11, Fillius et .al. C. 19 7 5 J observed a dramatic 
increase -- by a facto?: of thirty -- in the electron flux at 
E > 0.46 Mev, an energy compatible with the interpretation of the 
sheath model by Shawhan et al ■ Cl973, 19741. This "spike”, however, 
occurred outside the position of lo's orbit as inferred from the 
location of the sweeping minimum. This is in contrast to the 
prediction of the sheath model, in which electrons are accelerated 
away from Io by the Debye sheath on the inward face of Io and 
accelerated toward Io, striking its surface, in the outward face. 


Thus, it is not clear that the evidence for particle acceleration 
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at lo may be regarded as favorable evidence for the sheath model, 
although there exists the possibility that the combined effects 
of drift motions and diffusion may result in a displacement of the 
spike across the sweeping minimum. At this writing, however, the 
origin of these flux enhancements at Io seems ambiguous. 

III. THE INSTABILITY MECHANISM 


There is some variability in the literature in the nomencla- 
ture used for the various normal modes of plane waves in plasma. 
Therefore, we shall first define the nomenclature to be used here. 
We shall not be concerned with low-frequency modes for which the 
ion motion is important. In the high-frequency regime there are 
two normal modes in a uniform magnetized plasma. These are 
designated as the ordinary (o) and extraordinary (x) modes. The 

refractive indices n (= ck/w) of these modes are given by the 

o,x 


well-known Appleton-Hartree formula: 

2UJ 2 n - (JJ 2 /UJ 2 ) 

2 e v e' ' 

^ = 1 " I ^ 2 -* 2 ) Tfi 2 sin 2 0 ± Q A ’ 


(III-l) 


where UJ = 2ff , & = 2N £_ a , 6 is the angle between the wave 

e p,e e c,e’ ° 

vector k and the ambient magnetic field B q , and 

A — sin^ 6 + 4 (<JU 2 - tUg") cos 2 6j^. The upper sign refers to 
the o-mode, and the lower sign to the x-mode. If variations in the 
magnetic field or plasma density occur on a length scale much 
greater than k“^, Eq. (III-l) may be taken to describe the local 
propagation of the waves. 
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The ratio e = cu /fl is an important parameter in determining 
the nature of the normal -mode propagation. Throughout the region of 
the Jovian plasmasphere of interest for DAM,, e « 1. This is an 
unusual parameter regime in both laboratory and asirophysical plasmas; 
the only other known region in the solar system in which e « 1 is 
the high- latitude region above the terrestrial plasmasphere. Inter- 
estingly enough* this is the source region of the terrestrial kilo- 
metric radiation [ Gurnet t, 197^ j Kaiser and Stone; 19751* 

Tn Figure 8 we show the topology of the dispersion curves 

2 

uj versus k and the refractive indices n versus o> of the o- and x- 
modes for propagation at large angles to the magnetic field. The x-mode 
is seen to consist of two branches, the slow branch (over most of which 

9 n 

n x > 1) and the fast branch (n£ < 1). Both modes exhibit cutoffs, 

occurring at U) for the o-mode and u) = ^ {Ll + 2® ± (1 + 4®) 2 3/ 

' e x e 

2}^ for the x-mode, where the upper (lower) sign refers to the 
fast (slow) branch. For the fast x-branch, the cutoff frequency 
is approximately given by = ^ e (l 4- ®) . In addition, the slow 
x**branch exhibits a resonance at UJ = ^_(1 -f ®/2 sin^ 9). 

At 0 = 90° the o-mode is linearly polarized and the x-mode is 

■+ ■+ 

right elliptically polarized, with E j. B . As 0 0, the x~ and o- 

modes go over continuously to right and left circular polarizations, 
respectively. The portion a-b of the slow branch is the whistler 
regime. 

The x-mode cannot propagate in the "stop band" between the 
slow -branch resonance and the fast-branch cutoff. The axes of Figure 
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8 are normalized to SI ; for a fixed wave frequency, increasing c a/Sl 

S c 

generally corresponds to moving towards lower magnetic fields. 

Therefore, we see that waves on the slow x -branch may not propagate 
from the Jovian plasmasphere into interplanetary space. Strictly 
speaking, the stop hand exists only in the cold-plasma limit; thermal 
effects remove the resonance and lead to a non-zero minimum in the 
group velocity of the slow branch [Aubrey, Bitoun, and Graff, 1970] 

The finite temperature of the plasma also implies a collision fre- 
quency v , however, and waves with extremely slow group velocity V 
c g 

are attenuated by collisions over a scale length V /v . If the 

g c 

group velocity is sufficiently small, or the spatial region corre- 
sponding to the stop hand (the "stop zone") is sufficiently extensive, 
the attenuation in traversing the stop zone may still be large. Thus, 
as a general rule the existence of the stop band must be taken into 
account in assessing the propagation of waves. In particular, 
instability mechanisms in which the waves are amplified on the slow 
branch also imply the necessity of some mode-coupling process leading 
to the escape of waves in either the o-mode or the fast x-mode. 

The cold-plasma approximation may be deduced from the more 
rigorous kinetic theory in the asymptotic limit of low temperatures, 
under the conditions 

Other thermal effects on the dispersion relations are neglig- 
ible for our purposes, provided V^/c^ « e. Although e « 1, the low 
temperature of the Jovian plasma ensures the validity of this inequality. 
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(lll-2b) 


where k = k sin 0, k = k cos 0 . Equation ( III~2a) is the condition 

X II 

that finite-Larmor-radius effects may he neglected. Equation (ill -2b) 
is the condition that the thermal electrons not exhibit Landau and 
cyclotron resonance with the wave as they spiral along field lines. 

Such resonances may quickly damp the wave. 

Inequalities (III-2) arise from the inclusion of thermal effects 
in the plasma dielectric tensor. In the cold -plasma limit the 
dielectric is Hermitian, giving refractive indices which are either 
purely real (oscillatory solutions) or purely imaginary (exponentially 
damped solutions). Thermal effects modify somewhat the Hermitian 
part of the dielectric, which determines the normal-mode frequencies. 

As noted above, however, these modifications are of order V^/c^ for 
the electromagnetic modes, and are generality negligible. The thermal 
effects also introduce a non-Hermitian part to the dielectric. It 
is this part which implies damping of the normal modes; the damping is 
weak under the conditions (III-2). 

In the same manner, the inclusion of a population of energetic 
particles also modifies the plasma dielectric, in general introducing 
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a frequency shift into the real part of the refractive index and, 
in addition, an imaginary part to the frequency, which may then tie 
written 


ui(k) = aj r (k) + iy(k, u> r ) 


( III -3) 


The usual approach adopted in linear instability analysis is to 
consider an energetic particle population of low density relative to 
the thermal background plasma. Then the real frequency u) is governed 
by the cold-plasma dispersion relation,- while the growth or damping 
of the waves is determined by the energetic particle distribution. 

The dispersion relation is then of the form 


D(k, ud) “ C(l, c o) + £ <») 


(III-4) 


where C(k, to) is the cold-plasma dispersion relation and g s (k, u>) 
represents the contribution of energetic particles of species s. In 
principle, the thermal damping effects may be included in 2 s g g , but 
this Is not usually done; in particular, it is not done in the 
theories we shall discuss here. Therefore, the tendency for insta- 
bility must be compared a posteriori with conditions (lXI-2), to 
determine that the Instability is stronger than the competing damping 


by thermal particles. 


41 


All of the instability mechanisms we shall consider depend on 
an anisotropic distribution of a small energetic -particle population. 
Instabilities may exist in various parameter regimes in both w and 
in the parameters that describe the anisotropy , e.g. , T /T . The 
nature of the instability may be closely related to the nature of 
the coupling, and some coupling mechanisms are compatible with only 
certain types of instability mechanisms. To some extent, individual 
coupling mechanisms have generally been proposed with some particular 
type of consequent instability mechanism in mind, even when the 
coupling and instability have not both been investigated by the same 
author. Conversely, in some cases a proposed instability mechanism 
includes the tacit assumption that there exists an unspecified 
coupling mechanism which produces either the required anisotropy 
in the energetic-particle distribution or the conditions necessary 
for an existing marginally stable anisotropy to be driven unstable 
to electromagnetic wave amplification. 

It is sometimes asserted in the literature that the DAM 
emission mechanism must be cyclotron emission. There is no compelling 
evidence, either observational or theoretical, that this is the case. 
Goldreich and Lynden-Bell argue for coherent cyclotron emission pri- 
marily on the basis of the sharp beaming pattern of the emission, 
and also on the high intensity. Such beaming, however, Is an inherent 
feature of all of the mechanisms discussed below. In addition, there 
is the complicating factor that, if amplification occurs on the slow 
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x-branch, the beaming pattern may be closely related to or determined 
by the mode "coupling process by which the DAM escapes from Jupiter. 
Thus, cyclotron emission is only one of the possibilities for the 
instability mechanism. 

Before discussing individual instability mechanisms, it should 
be noted that a complete understanding of DAM may involve more than 
one mechanism. There are at least two distinct types of Io-modulated 
DAM, the decasecond (L) and millisecond (S) bursts. They have quite 
distinct time scales— from about 0.5-3 sec for L-bursts and 1-10 msec 
for S-bursts — which do not merge continuously into one another. 
Moreover, the S-bursts exhibit true frequency drifts, which are always 
negative. In contrast, L-bursts exhibit frequency drift only in their 
overall envelope ; this is probably a geometric effect of the beaming 
pattern. Although both L- and S-bursts share the same general 
morphology of Io -modulation, these facts imply that they may be caused 
by separate instability mechanisms. 

It is also unclear whether the Io-modulated and Io- independent 
L-bursts imply different emission mechanisms. With the recent recog- 
nition of an Io-independent component of the early source region 
LCarr and Desch, 1975; Bozyan and Douglas, 1975J, one might be tempted 
to hypothesize that the Io-modulated DAM is simply an enhancement of 
a weak continuous emission. It would then seem necessary, however, 
to shift from the common interpretation that the morphology of the 
I o -modulation implies a beaming pattern, and instead to suppose 
that the lo-moduiation is indeed only possible for certain 
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geometrical configurations relative to the observer. This latter 
hypothesis is physically untenable. Moreover, the declination effect 
in CML is related to the lo- independent DAM [.Goerts, 197l3, while 
the modulated component exhibits a declination effect in Io phase 
Lconseil, 1972; Lecacheux, 1974; Thieman, Smith, and May, 1975; 

Bozyan and Douglas, 1975 J. Thus, it seems probable that the lo- 
independent and -modulated components are due to distinct mechanisms, 
although these arguments do not conclusively rule out the possibility 
that they may be closely related. 

Thus, in assessing the viability of a given instability mech' 
anism, one should attempt to do so in the several contexts of Io- 
independent L~bursts and Io-modulated L- and S-bursts. In the 
first of these contexts, the lack of a demonstrable coupling mech“ 
anism is irrelevant . 

Finally, we note that the mechanisms we shall discuss have 
generally been analyzed only in the linear regime. Linear instability 
analysis may be quite misleading, however, because it provides no 
information on the saturation level of the instability. In general, 
a complete stability analysis must include consideration of the non- 
linear regime. 


A. The Sweeping Model 

Of the instability mechanisms we discuss, the sweeping model 
of Wu Ll973] is unique in that the instability is driven by aniso- 
tropy in configuration space rather than in velocity space. As was 
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mentioned in section II, the density gradients of trapped 

energetic particles across the boundary of the IFT produce drift 
currents in the cavity "wall 11 ; these currents are perpendicular to 
both B Q and N. Wu considered the motion of each energetic particle 

•"J 

species transverse to B Q to be characterized locally by a mean energy 
E ljS * The resulting drift velocity of the s-fch species is then given 
by 


<v 


> 

x,s 



(III-5) 


where q s is the charge, the density gradient is taken to be positive 
in the y-direction, and the drift motion is in the 4-x- direct ion for 
electrons and the -x-direction for protons. 

Wu considered instability at the upper hybrid frequency 
and for very short wavelengths, clc/^ £ >:> 1. Under these conditions, 
the energetic electrons contribute only a negligible shift in the 
real frequency, and the instability is produced by the non-resonant 
contribution of the protons. Taking an isotropic gaussian for the 
local distribution of the energetic protons, Wu found the growth 
rate of waves with k j_ B Q to be given by 
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where subscripts p,e denote energetic protons and thermal electrons, 

i- 

respectively, Rp = <E X 9 pM) z /^p , and is the x-component of the 
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wave-vector. Maximum growth occurs at 


k v_ 
x D 

fi 
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Taking , where 6 is the thickness of the density- 

gradient layer (see section II -C) , one finds v^ in the range 1-4 x 10 

km/sec for 5 MeV ^ < 10 MeV. Thus, the wavenumber at maximum 

growth is of order cK/^ ^ 100. Wu did not address the question of 

e 

how the waves penetrated the stop band. 

— » — ♦ 

Wu’s analysis was done for k strictly perpendicular to B q . 

Because the analysis was a local one, it also requires a finite 

component of k in the y-direction, so that 
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Although it is clear from Eq. (Ill -6) that Y increases as k /k 
increases, the inequality (III-7) is not sufficient to allow a precise 
determination of the azimuthal width of the beaming pattern; at 
best, one can say only chat it is a pencil beam. In the direction 
parallel to B q j the width of the pencil beam is sharply delimited 
by Eq. (III-2b) , which requires 
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The current -driven instability also contains the feature of 
a low-frequency cutoff, because an essential requirement for the 
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instability to occur is that, locally, e > 2m/M ^ 10 . In general, 

e decreases with increasing altitude, at least until the point of 
maximum zenopotential height [.Melrose, 1967; Smith, 1973; Goerfcz, 
19753. 

The power available in this mechanism was estimated by Wu 
[l973] to be about 10® W, adequate for the DAM requirement* Smith 
[1973] estimated the power to be 10® T|: W, where T| « 1 i s an effi- 
ciency factor. Smith also pointed out that the estimated require- 
ment of ~ lO^ 5 W depends on an assumed beaming into a solid angle of 
approximately 0.1 sr; from the above discussion this solid angle 
might be much too large in the context of the sweeping instability. 
Both of these estimates assumed growth rates of order Y 10 ~^ sec 1. 
Because the growth rate depends on the energetic-proton density, the 
precise value of the estimated power may have little significance. 

The ability of the mechanism to meet the power requirement of Io- 
modulated DAM, however, seems marginal. 

Comparison of the predictions and requirements of the sweeping 
model with Pioneer observations tends to cast doubt on the viability 
of the mechanism for Io-modulated DAM. As was noted above, sweeping 
seems to occur effectively enough to establish the necessary density 
gradient. Because the mechanism relies on trapped protons, however, 
and the higher-order moments of the magnetic field make the surface 

I 

field highest in the northern hemisphere, the instability is limited 
at any longitude to the region between the topside Jovian ionosphere 
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in the southern hemisphere and the conjugate mirror point. Because 
the growth rate of the instability is proportional to , the 
emission is essentially limited to the southern hemisphere. In order 
to span the range of DAM frequencies, therefore, the peak ionospheric 
electron density, N eo , must be sufficient to raise the upper hybrid 
frequency In the ionosphere up to about 40 MHz. The maximum gyro- 
frequency at the foot of the IFT is about 27 MHz, giving the require- 
ment that N e0 ^ 10 ^ cm - ^. The maximum density inferred by Fjeldbo 
et al . Ll975J, however, is of order 3 x 10~* cm ^ . Thus, unless this 
determination is too low by a factor of 30, the sweeping mechanism 
is essentially limited to frequencies below 30 MHz. 

Other aspects of the sweeping mechanism will be discussed 
below. In particular, it may be quite relevant to the question of 
the "lo-independent" DAM. 


B. "Maser -Like" Mechanisms 

A number of authors have considered wave amplification by 
electron distributions which are anisotropic in velocity space. Such 
mechansic are often called "maser-like”, because the anisotropic 
distribution is analogous to an inverted level population. The most 
commonly analyzed situation is that in which tut mea-' p st, endicular 
energy T l is greater than the mean parallel energy, Tjj . 

Chang Ll963] discussed the amplification of x-mode waves by 
an anisotropic distribution of relativisitic electrons. Although he 
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referred to these waves as whistlers, his analysis was in fact valid 
for the entire x-mode, He only considered the cases of propagation 
parallel and anti -parallel to 1^, however. Chang found two para- 
meter regimes in which amplification could occur. Of these, the 
only case relevant to the Jovian exosphere, in which e « 1, is 
the frequency regime 


2 2 
m » 0 m “ OJ > 0 



( III-8) 


Because € « 1, however, inequality (III-8) is equivalent to the 
condition 
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Thus, as in the case of the sweeping model, Chang's mechqnism requires 
ek/Q e » 1 (see Figure 8 ). For parallel propagation, however, these 
waves are strongly damped by the thermal electrons because, for 
^ « € and e « 1, we have 
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in violation of inequality (lll-2b). 
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In an analysis complementary to that of Chang, Goldstein, and 
Eviatar [1972] discussed amplification at 6 - 90 * on the slow x- 
hranch hy a loss-cone distribution of trapped electrons. Goldstein 
and Eviatar, however, explicitly recognized the necessity for the 
waves to penetrate the stop zone. They argued that the amplification 
would be sufficiently strong so that the radiation would simply 
tunnel through the stop band, suffering exponential attenuation but 
emerging into the outer region of the exosphere with observable 
intensity. Ihey asserted that this mechanism could be responsible 
for the lo- independent DAM. As we shall discuss below, the stop 
zone is probably sufficiently thin that exponential attenuation of 
waves incident upon it is negligible. Assuming this to be so, 
the instability analysis of Goldstein and Eviatar might account for 
lo- independent DAM by emission from a loss-cone distribution near 
marginal stability. 

Hirshfield and Bekefi L1963] considered the amplif ication 
of x-mode waves at ft , with 6 = 90°. In their model the distribu- 
tion function was assumed to be isotropic; the M maser-like M effects 
were due to the dependence of the gyrofrequency on the relativisitic 
electron mass. This dependence leads to coherent interaction 
between the particles and the wave field through phase synchroniza- 
tion, such as occurs in synchrotron accelerators [-Kuckes and 
Sudan, 197l3. 
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Although it was apparently not recognized by Hirshfield 
and Belcefi, the unstable waves at the gyrofrequency are below the 
stop band. They also predicted amplification at higher harmonics 
of the gyrofrequency, but at unobservable levels. Moreover, their 
work was done before the dis’cov'hry of the lo -modulation, and so the 
distribution function used in their analysis was based on the 
assumption of direct access to the -Jovian magnetosphere of streams 
of low-energy particles of solar origin. 

Goldreich and Lynden-Bell also postulated coherent cyclotron 
amplification for the emission mechanism, but recognized the problem 
of traversing the stop zone. They discussed the instability as seen 
in the frame of the streaming particles; in this frame their results 
are essentially similar to those of Hirshfield and Bekefi. Goldreich 
and Lynden-Bell asserted that in Jupiter's frame the waves would be 
Doppler-shifted above the stop band and so would escape the magneto- 
sphere. This conclusion was erroneous, however, because a wave 
with phase velocity smaller than c in one frame will not have phase 
velocity greater than c in any other frame. Thus, the Doppler-shifted 
waves postulated by Goldreich and Lynden-Bell would not be normal 
modes in the magnetospheric plasma, and would be heavily damped 
away from the source region. 

All of the mechanisms so far discussed provide for amplifica- 
tion on the slow -x-branch. With the exception of the sweeping -mod el 
instability, which is non-resonant and is carried by protons, 
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these mechanisms all invoke resonant-electron effects These are 
possible only on the slow x-branch, where the phase velocities of 
the waves may be exceeded by particle velocities. For the o-mode 
and the fast x-branch, however, the phase velocities are greater 
than c, and resonant interaction is impossible. There are, nonetheless, 
maser-like mechanisms to amplify these waves. ; they are sometimes called 
subluminous instabilities. Such mechanisms were hypothesized by Ellis 
[ 19 62 , 19 63 , 19 65 1 and by Ellis and McCulloch L 19631, who assumed the 
existence of electron "bunches'* precipitating into the Jovian ionosphere. 

The mechanism was investigated quantitatively in a series of papers 
by Fung [l966a,b,c]. The first of these papers [.Fung, 1966a] was 
concerned with terrestrial VLF emission and so neglected terms of order 

which are large in the Jovian exosphere. Fung 11966b] extended 
e e “ 

the analysis to include these terms. Both of these analyses assumed 
precipitating electron distributions of the form 

F (p) ~ 6(p - p ) 6(p - p ), (III-9) 

e I! ||o x xo 

and allowed for arbitrary angles of wave propagation. In his third paper, 
Fung L 19 66c] extended the formalism to include finite temperatures T^ 
and T|| in the precipitating stream. 

In this mechanism the emission is Doppler-shifted to a normal 
mode above the stop-band by the streaming velocity of the particles. 
Radiation at the fundamental is confined to a finite cone, in the forward 
direction, of wave-vector directions about the magnetic field, while 
emission at higher harmonics may occur in all directions but is much weaker. 
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Melrose [1973, 1975] alsg considered subluminous instability 
from a streaming anisotropic distribution, which he assumed to be a 
shifted b i-Maxwellian : 
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Here 0 = v/c, a is the particle pitch angle, and 0 = (0 cos a) is 

s 

the normalized streaming velocity. The Doppler- shifted frequency of 
the s-th hazmonic, radiated in the direction 0 by a particle of 
pitch angle a and Lorentz factor 7, is 


e 

7(1 - n (3 cos a. cos 0 ) 


At the fundamental, the condition that the radiated wave, be Doppler- 
shifted above the stop band is then 


n 0 cos 0 > e 
X s 


which is not excessively stringent in view of the smallness of e. 
Melrose finds a maximum growth rate of order 
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where T| = ^ Q /^ e is the ratio of the stream density to that of the 
background plasma. The bandwidth of the growing waves is of order 
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x 



« 1 
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which is in qualitative accord with observations. 

Melrose finds a threshold condition on the degree of aniso- 
tropy necessary in order for amplification to occur; the condition 
is 

£ 2 > 0 . (I3T-12) 

Although this inequality may seem stringent, it is actually fairly 
easy to meet under certain models for the acceleration of the parti- 
cles. For instance, the accelerated spectrum in the sheath model 
is essentially a delta-function in energy at any given point at Io; 
the particular injection energy depends on the spatial structure of 
the sheath potential [Smith and Goerta, 19751* Because the energy 
is in general much larger than the temperature of the particles before 
acceleration, the accelerated distribution will map at high latitudes 
essentially to a distribution of the form (III- 9 ) assumed by Fung. 
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The parallel motion may be equated locally to Melrose's p , while 
his may be equated to the characteristic transverse velocity. 

In the limit of this model, inequality (III-12) is satisfied virtually 
trivially. This condition appears to be an artifact of the assumed 
form (III-IO) for the distribution function; Fung [1966b] found 
instability for distributions of the -form (ill -9) in which < P * 
Melrose [19751 used qvasilinear arguments to estimate the 
saturation energy density W of the amplified waves as 





where is the gradient scale length of the magnetic field and 

-D 

2 2 2 

W = N me ( P + p )/2 is the energy density in the. stream. He 
S S J.O S' 

evaluated this for two different sets of beam parameters, obtaining 

-12 - 2 * * * * 7 

results which bracket an observational estimate of 4 x 10 erg cm 
required. This latter estimate is made using a very narrow beaming 
cone of solid angle 0.01 sr. In general, however, the energetic 
requirement of DAM appears to be well met within the context of 
Melrose's calculations. 

The maximum growth rates predicted by Melrose's theory are of 
2 3 -1 

the order of 10 ^ y < 10 sec . Thus the mechanism may be most 
relevant to L-bursts, rather than to S-bursts, which are already 
fully developed on this time scale. 
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Finally, we note that Melrose did not quantitatively address 
the question of the beaming pattern of his mechanism. As we shall 
discuss in section IV, this is an important topic for comparison of 
the predictions of the mechanism with the observational phenomenology. 

We have discussed only mechanisms that have been treated 
quantitatively in the literature. Other suggestions have been 
advanced by Warwick [196.1, 19633, Zheleznyahov [1966], and Marshall 
and Libby, For a discussion of these works, the reader is referred 
to the review by Warwick [1967]. 
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IV* WAVE PROPAGATION IN THE JOVIAN PLASMASPHERE 

Compared with the aspects of coupling and instability mech- 
anisms , propagation of the DAM waves in the Jovian plasmasphere has 
received relatively little attention. It is an important topic, 
however, which must he considered as an integral part of a complete 
theory of DAM. The beaming pattern implied by the DAM morphology 
is undoubtedly due to a combination of the inherent beaming pattern 
of the instability mechanism and the subsequent propagation of the 
waves. Conversely, the phenomenology may provide the basis for 
inferences regarding the plasma medium and the source regions. 

As an example of the latter approach, some authors have .on- 
sidered the polarization. Warwick and Dulk [1964] argued that the 
highly elliptical polarization of decametric waves upon emergence 
from the Jovian magnetosphere indicates that the base mode of the 
emission is the extraordinary mode 1 , 'ah ^interpretation which has 
generally been accepted in the literature. The sense of polariza- 
tion of this mode is right ( left) -handed when the wave vector has 
a component parallel ( anti parallel) to Although the limiting 

polarization depends upon the total ray path, the former case would 
generally be expected to occur for emission from the northern hemi- 
sphere, and the latter for emission from the southern hemisphere. 

At high frequencies (> 18 MHz), DAM is predominantly right-polarized, 
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while at lower frequencies the degree of left -polar! sat ion increases 
as the frequency decreases. The axial vector exhibits a quasi- 
sinusoidal rotation profile , which is conventionally interpreted as 
indicating that at a given longitude one or the other hemisphere 
dominates the beaming into a given direction. This reasoning should 
be considered with caution, however, for within storms that are rich 
in time structures there may exist considerable polarization diversity 
in the frequency-time plane. 

Notwithstanding the argument of Warwick and Dulk about the 
base mode of the emission, Warwick [1970] pointed out that the radia- 
tion incident upon the terrestrial ionosphere is not in a base mode. 
Nearly all of the Faraday rotation observed, however, is attributable 
to propagation through the terrestrial ionosphere; Warwick estimated 
that less than 360° of Faraday rotation occurs in the Jovian magneto- 
sphere. Warwick argued that this fact placed a small upper limit 
on the column density J N e ds through which the radiation propagates 
at Jupiter, and that this limit on the column density was incompat- 
ible with density models such as those of Gledhill [1967] and Goertz 
[1975, 1975], in which there is a disk-shaped concentration of plasma 

-ft 

in the equatorial plane. Goertz [197^1 challenged this argument; 

* 

Although the density models of Goertz and Gledhill are topo- 
logically similar, they are derived from different considerations and 
so differ by orders of magnitude in the actual densities predicted at 
a given radial distance in the plasma disk. 
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he showed that., if the radiation is generated entirely in a base mode, 
then during propagation through a spatially varying medium it will 
couple to the other mode just sufficiently to preserve its polariza- 
tion, and so Faraday rotation could not occur. 

As we have discussed in section III, several proposed insta- 
bility mechanisms actually amplify craves on the slow x -branch; such 
waves cannot escape directly to free space. The general topology of 
the stop zones has been discussed by Liemohn [1973 1 and by Gulkis 
and Carr [1966]; these latter authors suggested that the dimension 
of the stop zone might be longitudinally asymmetric, so that the DAM 
beaming pattern — as manifested in the source morphology- -might be 
governed by the location of the thinnest stop zones. They noted, 
however, that the mechanism by which waves' could penetrate the ^cop 
zones was unclear. 

Qya [ 197^1 attempted to account for the escape of waves 

through the stop band by considering mode conversion of extraordinary 

waves propagating inward towards the "planet- The conversion occurs 

by the reflection of the o-mode at the plasma level to = to- There 

are a number of difficulties with the geometry inherent in Oya's 

model. It suffices to note, however, that if the peak ionospheric 

5 

electron density is 3 X 10 as reported by Fjeldbo ej^ al. , then the 
corresponding plasma frequency is about 5 MHz, and so the critical 
condition is not meet at DAM frequencies. 

The importance of the stop zones, however, is greatly depen- 
dent on the temperature and peak density of the Jovian ionosphere. 
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Smith. [1975] showed that for N o = 3 x 10^ [Fjeldbo et al. , 1975] 
and T q = 750° [KLiore et al. , 1975 lj the attenuation of waves during 
traversal of the evansecent region is negligible for frequencies 
helow about 20 MHz. This fact has profound consequences for the 
explanation of Io -independent DAM, which we discuss in section V. 

As we have seen in section III, theories in which instability 

occurs on the slow x-branch, such as those of Goldstein and Eviatar 

[1972] and of Wu [1973 1* amplify waves of very low group velocity 

V « c. Smith et al. [1972] and Wu et al. [19731 attempted to 

s 

account for the lo- phase asymmetry (referred to in section II) in 
terms of the low group velocity of the waves before their escape 
from the Jovian magnetosphere. In addition, these authors suggested 
that the sense of frequency drift of the modulation lanes — negative 
in the main and third sources and positive in the early source — was 
caused by he vectorial combination of a low group velocity in the 
medium relative to the observer. This latter suggestion stemmed 
from an interpretation of the modulation lanes as drifts as a modu- 
lated structure in the emission itself. Smith [19731 showed that 
this interpretation could not be correct if the modulation structure 
in the emission were impressed upon the IFT; instead, a longitudinal 
modulation is required within the framework of the group-delay 
hypothesis. 

Smith [1973] also found that refractive effects upon upper- 
hybrid waves were strong, tending to make the waves propagate into 
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the stop zone even when the wave vector is initially perpendicular 
to the no mal to the stop zone surface. In view of Smith's later 
results concerning the thickness of the stop zones, it is unlikely 
that group-delay effects play any significant role in either the 
phenomenology or morphology of DAM. 

One caveat to be made about this last statement is that it 
implicitly assumes geometrical optics to be valid everyvrhere . All 
of the instability mechanisms discussed in section III imply either 
subsequent propagation through a thin evanescent region or, as in the 
case of Melrose's mechanism, that the radiation is intially beamed 
toward the planet and subsequently reflected toward the observer. 

All studies of wave propagation to date have been in the context of 
geometrical optics, in which the medium is slowly varying and t’.e 
propagation at any point is governed by the Applet on -Hartree equation 
with the local plasma parameters. The condition for the validity of 
geometrical optics, which is essentially a WKB approach, is 

— ~ | Vn j « 1 . (IV-1) 

con 


Wear resonance and reflection, inequality ( IV-l) breaks down* Never- 
theless, application of geometrical optics in the ease of reflection, 
for example, usually yields good results for the ray trajectory even 
when the solution has been carried through a thin region where WKB 
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analysis Is invalid [Ginzburg, 19T0}. Strictly speaking, however, 
the breakdown of inequality ( IV-l) in some region implies that a 
full-wave solution must be carried out in that region and matched to 
the WKB solution in regions where geometrical optics is valid. It 
is possible that such an analysis would explain some of the fine 
structure observe !n decametric spectra. 
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V. SUMMARY AND DISCUSSION 

For the most part the discussion above has treated coupling, 
instability, and propagation separately, although occasional refer- 
ences have been made to the relations between them- In this section 
we attempt to draw the different aspects together in evaluating the 
present state of the theory. We shall also indicate some of the 
major unsolved problems, and suggest some possible approaches for 
future work. 


A. Io -Independent: DA M 

As we have already remarked, there are not one but three 
distinct decametric emissions: lo-independent L-bursts and Io-modulated 

L- and S-bursts. Since the discovery of the lo-modulation only one 
paper, that of Goldstein and Eviatar, has been addressed specifically 
to the question of an emission mechanism that cannot be related to a 
coupling mechanism. This mechanism appears to be viable for the Io~ 
independent DAM, especially in light of the indications that tunneling 
of the radiation through the evanescent regions does not lead to severe 
attenuation. The radiation may emanate from particles trapped on L- 
s hells close to but inside the orbit of lo, and so may be unaffected by 
any magnetospheric perturbations local to Io. The waves are inherently 
beamed sharply perpendicular to the magnetic field, in accord with the 
declination effect. The chief shortcoming of the work of Goldstein and 
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Eviatar is that they did not estimate the level at which the instability 
would saturate due to pitch-angle scattering by the radiation,, 

Another theory which may account for the Io- independent DAM 
is the sweeping model of Wu, Although this model was originally 
proposed to explain lo-modulated emission, there are several reasons 
why it may be more relevant to the unmodulated radiation. Among these 
are the sharp beaming pattern perpendicular to the magnetic field, 
and the fact that, as discussed in section III, the instability is 
limited to frequencies below 30 MHz within the context of the density 
inferred by Fjeldbe et al . In addition, the growth rates evaluated by 
Wu and by Smith [l973] are rather low, and cannot account for the time 
scales of lo-modulated L-bursts. 

Although sweeping is certainly a coupling mechanism, the real 
question is whether its effects are localized close to Io. Smith C. 197 3 3 
showed that in addition to the upper-hybrid instability, the drift 
currents induced by sweeping also drive a long-wavelength instability 
at the lower hybrid frequency = + e(M/m.)/( 1 + e) . Smith 

evaluated the spatial diffusion coefficient [.Dupree, 1967] 


D 

JL 


max 
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and showed that it depended on both the thermal plasma density and 
the local gyrofrequency , in addition to its dependence upon the 
energetic-proton distribution. In general, the diffusion is slower 
at high altitudes, so that the density gradient is maintained longer 
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at low frequencies. Smith suggested that the lo- independent DAM is 
in fact due to sweeping, but is poorly correlated to the position of 
Io because of the slow diffusion at low frequencies, and thus appears 
statistically independent of Io. Another diffusion mechanism has 
recently been investigated by Huba and Wu C 19 75 ] . 

Wu Cl9 733 and Smith C 19733 neglected the longitudinal drift 
motion of the protons in the dynamics of the sweeping process. This 
drift is oppositely directed to lo's orbital motion relative to the 
corotating plasma. As shown by Thomsen and Goertz Cl975i, the swept 
cavity drifts in the same way. Thus, the effective sweeping time is 
reduced from 13 hours to about 4 hours, and it is likely that the cavity 
persists along the entire Io ortit. Although the diffusion tends to 
reduce the magnitude of the density gradient during this period, t-he 
gradient is never completely eliminated, and so the drif t-current-driven 
instability is always occurring. 

A possible impediment to the sweeping model as the cause of 
unmodulated DAM might exist if in fact Io causes local heating of the 
ionosphere by MHD waves, as suggested in section II. Such heating 
would drive thermal plasma up into the IFT, increasing the scale height 
H of the thermal -plasma density gradient. If H» eEjj then the thickness 
of the stop zone becomes equal to H, so that x-mode waves could not 
penetrate the stop zone. We might expect this enhanced density '’filament" 
to cool quickly by thermal conduction to the ionosphere, however, and 


so be localized close to Io. 
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B. lo-Modulated L-Bursts 

A complete theory for lo-modulated DAM would synthesize all 
three aspects of the problem,, Starting from a particular demonstrated 
coupling mechanism, the theory would predict the form of the resulting 
free-energy source, for example the velocity-space distribution of 
accelerated particles. It would then incorporate this distribution 
into an instability analysis, demonstrating that the instability could 
account for the various aspects of DAM phenomenology such as time scales, 
bandwidth, nolarization, and spectral power. The inherent beaming pattern 
of the radiation in the source region would be evaluated, and its 
subsequent propagation in and emergence from the Jovian plasmasphere 
would be compared with DAM morphology. In addition, the theory should 
contain some natural explanation for the probabilistic nature of the 
emission. 

Needless to say, such a comprehensive theory has not yet been 
attained. Nevertheless, some promising work has been done on different 
aspects of the problem. Smith and Goertz have given the sheath model 
of acceleration a firm theoretical basis. The emission mechanism of 
Melrose is compatible with the sheath coupling model, and seems able 
to account for most of the necessary features of the phenomenology. 

It yet remains, however, to synthesize these two models, and to 
investigate the propagation and reflection of the resulting radiation. 

Continued investigations of other mechanisms for both coupling 
and instability might also yield viable new models. For example, 
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Palmadesso ej: al. C 197 5 3 have recently proposed a parametric insta- 
bility mechanism for the terrestrial lcilometric radiation. This 
theory seems attractive, and might also have application in the 
parameter regime of DAM. 

There can be no qtte-stion that propagation effects contribute 
significantly to the DAM morphology, and pttfbably also to the 
phenomenology. Any instability mefchanism has an inherent beaming 
pattern in the source region. The structures of the magnetic field 
and the exospheric density distribution combine to modify this pat- 
tern into the directivity pattern that is conventionally inferred 
from the morphology. A complete theory of DAM must include a study 
of the "antenna pattern" of the Jovian plasmasphere to show that it 
is compatible with this morphology. Such a study would serve two 
purposes. First, to the extent that it is dependent upon a model of 
the exospheric structure, it would help to elucidate that structure. 
Second, the emerging radiation pattern depends upon the beaming pattern 
of the instability in the source region. Thus, study of the propaga- 
tion may place boundary conditions upon the instability mechanism 
which allow one to discriminate between likely mechanisms. Because 
the instability mechanism is, in general, also dependent upon 
the ambient plasma regime, we might expect that the synthetic study 
of coupling, instability, and wave propagation will lead to a con- 
vergent picture. 
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C. Three Major Problems 

None of the models we have discussed appear able to account 
for the lo-modulated S-bursts. Ellis Cl965] interpreted them as 
being due to streams of 30-60 keV electrons moving outward. It is 
unclear why such streams should always be accelerated outward, or 
why they should occur primarily in the early-source region. It is 
also unclear why they should give rise to a different instability 
than that responsible for the L-bursts° As noted in section III, 
the well-separated time scales of the L- and S-bursts, and the time 
frequency drifts of the S-bursts, indicate that different mechanisms 
produce the decasecond and millisecond emissions. For the moment 
the S-bursts remain enigmatic. 

A second major question is the explanation of the asymmetry 
in Io phase for emission viewed from east and west of the Earth-Jupiter 
line. This may be a local-time propagation effect; if so, its explana- 
tion would have important implications about the structure of the 
Jovian plasmasphere. 

Finally, one of the most intriguing problems regarding DAM 
is the explanation of the modulation lanes observed in both L- and 
S-bursts. Riihimaa C 19 71, 19743 has documented an extensive morphology 
for these features. The modulation lanes probably contain a great 
deal of information about the mechansim and propagation of DAM. 



Epilogue 


We began this review with the observation that a comprehensive 
understanding of DAM has yet to be attained. Nevertheless, many 
authors have contributed a great deal of insight to our perspective 
of the problem. The conceptual framework for the theory is largely 
in place, and we may reasonably expect much progress in the future. 
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APPENDIX A 

HEATING OF THE JOVIAN IONOSPHERE BY MHD WAVES 


In section II-A we mentioned the possibility that MHD waves 
generated by Io might provide local heating of the Jovian ionosphere. 
This comment was based on the following rough but indicative calcu- 
lation, due to Wu and Smith [unpublished]. 

Assume the waves steepen into shocks during their propagation 
along the IFT. The typical dimension of a eollis ionless shock is 
c/m. , where ci>. = Sn f . is the angular ion plasma frequency. Then 

' i i pi 

for a wave amplitude of SB, the current density j is given by 


m, SB 

5 - 5S h X B ! - 


The Joule heating rate of the plasma is then 



.2 

cr 


2,^2 

a>. SB 

l 

l6h 2 CT 


(A-l) 


The conductivity a can be either anomalous or the usual Coulomb 
conductivity 
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KgT 2 

CT c 2 » 

m s e te A 
e 

where A = (K^T) ? / (4 ttIT) s e . Use of a Q will lead to an upper 
hound for the heating. Substituting (A-2) with (A-l), we obtain 


(A-2) 


B_ t 5/2 _ 5_ 
at 32 


A 2 2 

m e s e Bn A tu^ 6B 

Vh*/'*'"' 


(A-3) 


Because the temperature dependence on the RHS (A-3) is logarithmic , 
we may obtain an approximate solution by taking faA to be constant j 
we assume Qm.A^ 20. Assuming that all the heating occurs in the 
ionosphere, in a time interval T , we find 


T(r ) - T q 



13 

1.7 x IQ ^ 
T 5 / 2 



(A-4) 


Taking T q =- 750 K [Kliore et al. , 19751 and. r ^ 60 s (the time 
for which the IFT threads lo), (A-4) becomes 


T(t = 60 sec) *- 


10 6 6B^ 5 


(A-5) 


We may estimate 6B in terms of the amplitude 5B^ of the 
perturbation by Io. For a rough estimate, we consider a generalized 
Pointing theorem. Let U be the total energy density (electric, 
magnetic, and kinetic) in the wave, and let ^ denote the group 
velocity. The generalized Poynting theorem is: 
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^ + V - VU + UV-V = 0 . (a-6) 

St g g 


We assume that a WKB solution is possible. (Strictly speaking, this 
is not consistent with our estimate of |V x "S| from the shock condi- 
tion; we discuss this below. ) The Alfven wave dispersion relation is 


(JO 



therefore , 


V = ctB , 


2 2 _ i 

where cl = [1+rrWm (l + V./c )] e , and so 
P A 


V * V = B(s) ^ a(s) = f(s) 


— f 

where 3a/cis is the directional derivative along B. Assuming a solu- 
tion of the form 


U - W exp - J 


„ f(x) 

vTST 

£ 


(A-7) 


Eq. (A-6) transforms to 


the solution of which is simply 


W( s, t) = W 




The boundary condition on the solution is 


w(o, t) = u Io (t) - 


Furthermore , the integrand in (A-7) is singly 


f(s) 


a 

3s 


i 


so that 


a. 


U < s > ‘ u Ioa^>- u Io 


N 


Io 1 + 7 


(A-9) 


In an Alfven wave, energy is partitioned in the ratios 
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E B 1 tt 2 
8tt : 8tt : 2 pV 



1 


2 2 

Thus, neglecting V^/c throughout, we may combine { A— 5) and (A-9) 
to find 


T(t- ) 


10 


'M 


2/5 


(6B ro J 4/5 




where N is the ionospheric electron density- Assuming W =- 30 cm 
o Xo 

and N q — 3 X 10^ cm’"" 5 , we have T(t) =« 4 X 10^ (fiB The power 

radiated by lo in Alfven waves is approximately 


V A,I 0 


f e 4 

\ 8tt 


Taking =* 7 X 10" 5 km/sec, we require fiBj Q =* 2 x lO’"' 5 g if 

10 -5 14 

P = 10 W, and 2 X 10 g if P = 10 W- (The former figure is 

Schmahl 's estimate of the power; the latter is the upper limit for 

which the effect on Io T s orbital period would not be observable.) 

3 

The corresponding temperature, according to (A-10) are 7 x 10^ K 
and 3 X 10^ K. 

The estimate (A-10) is quite crude, however, because the 
arguments leading to it were not entirely consistent. On the one 
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hand, we assumed that the waves steepen into shocks in order to esti- 
mate the Joule heating current* On the other hand, we assumed that 
the wave propagation was compatible with a WKB solution in order to 
find the steepened amplitude. The effect of the first assumption is 
to overestimate the numerical coefficient in Eg. (A-10), the effect 
of the second is to underestimate 6B. The value of the estimate 
should not be taken too seriously; rather, the calculation should be 
viewed in the spirit that the large discrepancy between the tempera- 
ture as given by (A-10) and the nominal ionospheric temperature of 
750 K indicates that transient local heating by lo may be significant. 
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FIGURE CAPTIONS 


Figure 1- The density of protons in the energy range 108 eV— 4.80 
keV near Io, as measured by Frank et al. [1975]* These 
authors fit their proton density measurements throughout the 
inner plasmasphere to Maxwellien distributions of character- 
istic temperature =* 100 Ey and densities between 50 and 100 
-3 

cm. 

Figure 2. The geometrical configuration of the Earth-Jupiter-Io 

system for peak occurrence probability in the various sources. 
Note the asymmetry in Io phase with respect to the Earth- 
Jupiter line. The CML ranges of the sources are in System 
III ( 1957*0) coordinates, at the Pioneer 11 epoch (1974-9)* 
Figure 3* (After Goertz and Deift [1973]-) 

(a) The magnetic field topology near Io. 

(b) Detail of the x-type neutral point reconnection region 
behind Io. 

Figure 4, Schematic of the Io sheath model [after Shawhan et_ al. 
[1974]). 

Figure 5- The Idealized circuit analyzed by Gurnett [1972]. The 
inset shows the plasma flow velocity in the context of 
Gurnett T s analysis. 2^ is the height-integrated Pedersen 
conductivity of the Jovian Ionosphere. For 2^ less than a 
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critical value £ ^ 5 mho, the electric field in the iono- 
c 

sphere is large and the flux tube tends to core ate with Io. 

For S > E , the voltage drop across the ionosphere is small 
P c 

and the plasma slips by Io. 

Figure 6 . Cross-section of the energetic-particle cavity in Io's 
orbital plane. 

Figure 7. Mapping of the density-gradient layer along the Io flux 

tube. Also shown are the directions of the density gradients 
and drift currents at the inner and cuter walls of the cavity. 
Figure 8 . Dispersion curves and indices of refraction of the ordi- 
nary and extraordinary modes for e « 1. The approximate 
regions of instability found by various authors are indicated. 
Figure 9* The gyrofrequency at the feet of the Io flux tube, in the 
magnetic field model of Smith et_ al. [1975 ]• The approximate 
sub-Io longitude ranges of tbe various sources are indicated 
(epoch 197*4-. 9 )* 



PROTONS (CM) 


i 


i 


ro 



Figure 1 

*** m iajm 



0° SGC 

i 

i 

i 

i 



0° SGC 




Figure 2 






COLLECTED 
THERMAL x 
ELECTRON FLUX 


JOVIAN MAGNETIC 
FIELD LINES 



10 IONOSPHERE 
CONDUCTIVITY 


POSITIVE 

SHEATH 


PLASMA 

POTENTIAL 


F igure 4 











z 



Figure 5 



DIFFUSION 


Figure 6 



2.5 



Figure 8(a) 



MELROSE 



v 

CHANG ( Q =0°) 

WU ( 6 =90°) 





f c (MHz) 






